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^ ; ABSTRACT 

We present the results of a targeted 3-mm spectral line survey towar ds the e i ghty-t hree 
6.67 GHz methanol maser selected star forming clumps observed by iPurcell (l2006h . In 
addition to the previously reported measurements of HCO+(1-0), H^^CO"'" (1 -0), and 
CH3CN (5-4) & (6-5), we used the Mopra antenna to detect emission Unes of N2H+ (1-0), 
HCN (1 -0) and HNC (1-0) towards 82/83 clumps (99 per cent), and CH3OH (2-1) towards 
^ ; 78/83 clumps (94 per cent). 

The molecular line data have been used to derive virial and LTE masses, rotational tem- 
peratures and chemical abundances in the clumps, and these properties have been compared 
between sub-samples associated with different indicators of evolution. The greatest differ- 
ences are found between clumps associated with 8.6 GHz radio emission, indicating the pres- 
I ence of an Ultra-Compact Hii region, and 'isolated' masers (without associated radio emis- 

sion), and between clumps exhibiting CH3CN emission and those without. In particular, ther- 
mal CH3OH is found to be brighter and more abundant in Ultra-Compact Hn (UCHii) regions 
k>( \ and in sources with detected CH3CN, and may constitute a crude molecular clock in single 

i dish observations. 

' Clumps associated with 8.6 GHz radio emission tend to be more massive and more lu- 

minous than clumps without radio emission. This is likely because the most massive clumps 
evolve so rapidly that a Huper-CompactHii or UCHii region is the first visible tracer of star- 
formation. 

The gas-mass to sub-mm/IR luminosity relation for the combined sample was found to 
be Lex M°-^^, considerably shallower than expected for massive main-sequence stars. This 
implies that the mass of the clumps is comparable to, or greater than the mass of the stellar 
content. 

We find also that the mass of the hot core is correlated with the mass of the clump in 
which it is embedded. 

Key words: ISM:molecules — stars: formation — ISM:abundances — surveys — stars:pre- 
main-sequence 



1 INTRODUCTION 

Molecular emission is a powerful tool when used to investigate the 
physical and chemical conditions in hot cores. Transitions requir- 
ing different temperatures and densities for excitation constitute an 
excellent probe of physical structure. Because the chemical proper- 



ties of hot cores vary with time, the relative molecular abundances 
can also be used as indicators of evolution. 

In the preceding decades, representative line surveys of a 
limited sample of cores have begun to identify the molecules 
most suited to investigating the process of massive star forma- 
tion. To date these have been restricted to a few objects: Orion- 
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G323.74-0.26 



Sutton et al.'l991^■ G34.3+0.15 jMacDonald et alJl996l : lKim etM 

2000, . 2001), IRAS 17470-2853 (also known as G5. 89-0.39, 



Kim et al.ll2002h . G305.2+0.2, jWalsh & Burtor]l2006h . 

Chemical models, including both gas-phase and grain-surface 
reactions have been developed for specific objects whose phys- 
ical structure and conditions are well k nown, e.g. G34.3+0.15 
jMillar et alj|l997l: [Thompson et alj|l999h , as well as for general 



Rodgers & Chamlev 1 2003h . These and earlier models 



cases, e.g. 

are limited by the computational power available at the time and 
consider either time-dependent chemistry at a test position, or a 
'snapshot' of spatial abundances at a single time. Advances in com- 
putational power have recently allowed the development of models 
which consider both time- and space-dep endent chemistry s imulta- 
neously, e.g., the model of AFGL 2591 bv lPotv et an i l2002h . Today 
it is the lack of observational constraints which restrict our under- 
standing of chemistry in young massive stellar objects. The priority 
in the coming years must be to assemble a large sample of molecu- 
lar abundance measurements towards a broad range of massive star 
forming regions, at different stages of evolution and with luminosi- 
ties and masses from B3 to 05 type massive stars. 

In this paper we present the remaining results of the Mopra 
'Hot Molecular Cores' (HMC) survey towards 6.67 GHz methanol 
(CH3OH ) maser selected targets. The initial results were pre- 
sented in IPurcelll ( [2003) (hereafter Paper 1), wherein we reported 
the detection of the hot-core species methyl-cyanide (CH3CN) 
and the observations of the formyl ion (HCO'*'), and its iso- 
topomer (H^''CO^). Here we expand the analysis to cover ther- 
mal transitions of methanol (CH3OH), carbon-monoxide (^"'CO), 
diazenylium (N2H"''), hydrocyanic acid (HCN) and hydroisocyanic 
acid (HNC). We analyse the full complement of molecules for evi- 
dence of evolution. 

Note: Tables 3-9 are available in electronic form only along 
with additional supporting text and figures. 
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Figure 1. Sample spectra for the maser site associated with G323.74— 0.26. 
The vertical dashed line marlcs the systemic velocity of the source, mea- 
sured from the peak of the optically thin H^'^CO^ line profile (see Pur cell 
i2Q06.) . Vertical dotted lines mark the expected velocity of other transitions 
of the same molecule. Individual Gaussian fits to the line profiles are plot- 
ted by a dotted grey line and the cumulative fit by a solid grey line. We did 
not attempt to fit Gaussians to HCN lines as the majority of profiles have 
complex shapes. 



2 OBSERVATIONS 

Observations were conducted on the Mopra millimetre wave tele- 
scope over five years 2000 - 2004 during the winter observing sea- 
son, lune - September. We observed the transitions CH3OH (2- 1), 
"CO(l-O), N2H+(l-0), HCN (1-0) and HNC (1-0) as single 
pointings, targeted at the maser sites. Electronic properties of the 
molecules and transitions are presented in Table [T] Five sources 
do not contain methanol masers but were detected as discrete 
clumps at 1.2 mm , adjacent to clumps containing CH3OH masers 
( iHill et al.ll2005h . These 'mm-only' sources were observed as po- 
tential precursors to the hot-core phase of massive starformation. In 
each case we tar geted the peak of th e 450 [im continuum emission 
later detected by IWalsh et alj 12003). The signal from the receiver 
was processed in an autocorrelator backend configured to have a 
bandwidth of 64 MHz, split into 1024 channels. At 3-mm wave- 
lengths this delivers a velocity resolution of ~0.2kms^^ over a 
usable velocity interval of ~130kms~^. 

All observations were performed in position switching mode 
by integrating for equal times on the science object and an 
emission-free reference position. Initial reduction of the raw data 
was performed using the SPcf^ package. Polynomial baselines were 

^ SPC is a spectral line reduction package for single dish telescopes, writ- 
ten by the ATNF. http://www.atnf.csiro.au/software/| 



subtracted from the line-free channels using the software pack- 
age. The CLAS^ spectral line analysis software was used for all 
further processing and analysis, including baseline subtraction and 
fitting Gaussian line profiles. For further details on the observing 
setup and data analysis please refer to Section 3 of Paper 1. 



3 RESULTS 

Thermal molecular emission was detected in 82/83 targets, with the 
exception being the source GIO.IO— 0.72. In the interim period this 
source has been reobserved at 6.67 GHz using the 30 m Ceduna 
radio telescope (University of Tasmania) and no CH3OH maser 
was detected (S. EUingsen, private com munication). This s ource 
is thus most likely a noise artifact in the IWalsh et alj ( Il998l) sur- 
vey and is discounted in the remainder of this work. Emission 
from ^^CO(l -0), N2H+ (1-0), HNC (1-0) and HCN (1 -0) was 
detected towards the remaining 82 sources, while CH30H(2-1) 
emission was detected towards 79/82 sources (96 per cent). The 

^ XS is written by P. Bergman, Onsala Space Observatory. 
ftp://yggdrasil.oso.chalmers.se/pub/xs/ 

CLASS is part of the GILDAS (Grenoble Image and Line Data Analysis 
Software) working group software. http://www.iram.fr/IRAMFR/GILDAS/| 
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3-(T noise limit on thie spectra where no lines were detected is 
recorded in Table[2] Figure[T]presents sample spectra for the source 
G323.74— 0.26. The vertical dashed line marks the systemic veloc- 
ity of the gas, measured from the optically thin H^^CO^ line (Pa- 
per 1). In the case of CH3OH, the rest-frame is centred on the A-l- 
transition, and for N2H"'" and HCN the rest-frame is centred on the 
Jf = I2 Oi E transition. Plots of all molecular spectra for 
all 83 observed sources (including GIO. 10—0.72) are available as 
additional online material. 



3.1 Line profile parameters 

As with CH3CN in Paper 1, we fit the multiple transitions of 
CH30H(2-1) simultaneously with four Gaussian lines having 
equal full-width half-maximum (FWHM) linewidths and whose 
line-centre separations were fixed to the theoretical values. The 
weakest component in the spectrum, the Jk.k' = 2i,i li,o E 
line, was detected towards 50 per cent of sources. Individual line 
profiles are generally well fit by single Gaussians, and the av- 
erage linewidth is 4.7kms~^. The CH3OH line profiles of five 
sources (GO.21+0.00, G0.26-h0.01, G6.54-0.11, G10.30-0.15, 
and GI5.03— 0.71) were too confused or too weak to be success- 
fully fit with Gaussians. The parameters of the Gaussian fits to the 
remaining 74 spectra are presented in Table 3. 

The N2H'*' (1 -0) transition is split i nto seven hyperfine lines, 
as illustrated in Figure |2] (adapted from lCaselh et alj|l995h . The 
relatively broad linewidths (> 2 km s~^) observed towards massive 
star forming regions result in the seven components blending into 
three groups with roughly Gaussian shapes. Consequently, we have 
characterised the profiles using two analysis schemes. In the first 
instance we used the hyperfine structure fitting (HFS) routine in 
the CLASS software to simultaneously fit the N2H^ profiles with 
seven Gaussians. Free parameters were the velocity (Vlsr), com- 
mon linewidth, excitation temperature, and optical depth. In the 
case of severely blended lines (as presented in this work) the fits are 
not well constrained and values for the optical depth and excitation 
temperature contain substantial uncertainties. Using this method 
we find the mean N2H"*' linewidth of the sample is 3.0 kms^^. 
Because of the uncertainty introduced by the blending, CLASS as- 
signed an effoneous value of r — 0.1 to all sources, rendering 
the resultant excitation temperatures meaningless. As an alterna- 
tive measure of optical depth and integrated intensity, we also fit 
the three blended groups individually with single Gaussians. In to- 
tal we successfully fit Gaussians to 79 spectra. The N2H"'" emis- 
sion in the source G0.26-l-0.01, G6.54-0.il and G6.61-0.8 was 
too weak or two confused to achieve a reliable fit. We discuss the 
optical depth and column densities derived from the N2H^ spectra 
in Section l4!4l The results of the CLASS HFS fits and the parame- 
ters of the Gaussian fits are presented in Table 4. 

Line profiles of HNC are observed to be similar to HCO^. 
The critical densities of the two molecules jire approximately the 
same at ~ 3 x 10^ cm"^ dSchilke et al.lll992t) and it is likely that 
they trace the same gas. Both molecules exhibit high velocity line- 
wings. Self-absorbed HNC line profiles are common in the sample 
and we referred to the low optical depth H^^CO''" and N2H+ line 
profiles when attempting to distinguish between self-absorption 
and multiple clouds along the line of sight. We fit self-absorbed 
profiles with a single Gaussian by masking off the absorption dip 
(see Paper 1, Figure 3). We also measured the line intensity by in- 
tegrating under the line between velocity limits bracketing the pos- 
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Figure 2. Spectrum of N2 H"*" (1-0) towards a dark cloud showing the 
seven hyperfine components jCasellietal.ll995l) . Linewidths towards mas- 
sive star-forming regions are 4-5 times greater, resulting in the hyper- 
fine components blending into three groups, designated 1 - 3 in this work 
(dashed lines). 



itive line emission and we calculated an equivalent linewidth from 



AVo 



/ Tmb dv 



(1) 



In Equation [T] J Tmb dv is the integrated intensity and Tmb is the 
peak main-beam brightness temperature. HNC emission from the 
sources GO. 21-1-0. 00 and G0.26-f0.01 was not measured as the line 
profiles were too confused. In total 79 spectra were successfully fit 
with Gaussians, while 80 were measured by integrating under the 
line profile. The results of the fits and integrated intensity measure- 
ments are presented in Table 5. 

The ^^CO line profiles often exhibit multiple components 
spread over the 64 MHz bandpass. Usually a single bright com- 
ponent is present at the systemic velocity of the source, which 
is easily fit with a single Gaussian. Occasionally several smaller 
components are crowded around the main line and these are fit as 
blended Gaussian lines. In general, the ^''CO profiles appear mod- 
erately optically thick but exhibit less self-absorption than HNC or 
HCO^. We again measure the integrated intensity and peak bright- 
ness temperature directly from the raw spectrum. The parameters 
of the fits and integrated intensity measurement of all 82 spectra 
are presented in Table 6. 

HCN exhibited highly confused line profiles, with few sources 
exhibiting the expected 3-component hyperfine structure. We at- 
tribute the complex profiles to a combination of high optical depths, 
multiple clouds along the line of sight and low level emission from 
broad line-wings tracing outflows. It is possible that some spectra 
may be contaminated by emission from the reference position, as 
HCN appears to be ubiquitous near the Galactic plane. 

Like HCO^, high velocity outflows wings are common in the 
spectra of HNC and, to a lesser extent, ^"^CO. We fit broad Gaus- 
sians to these wings simultaneously with the main component, and 
their parameters are presented in Tables 5 and 6 marked with a 'w'. 
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Figure 3. Distributions of full-width half-maximum linewidths from the Gaussian fits to the spectra. The HNC and ^''CO fine profiles were fit with single 
Gaussian components. CH3OH was fit with multiple Gaussians, where the linewidths were constrained to have the same value. The N2H+ spectrum was fit 
using the HFS method in CLASS. Seven Gaussian components with a common linewidth were fit to the three groups of blended lines. The value plotted here is 
the linewidth of an individual component. The mean CH3OH, HNC and ^^CO linewidth are all at approximately 4.5 km s~^, while the mean N2H+ linewidth 
is at 3.0kms^^. 



The Gaussian parameters for blended lines are marked with a 'b' in 
the same tables. 



3.2 Linewidths 

Figure |3] shows the distributions of full-width half-maximum 
linewidths for the four molecules which were fit with Gaussians. 
The median CH3OH, HNC and ^^CO linewidths are all at ap- 
proximately 4.5kms~^, while the median N2H+ linewidth is 
uniformly lower at 3.0 km s^^. Kolmogorov-Smirnov (KS) tests 
yield probabilities above 50 per cent that the CH3OH, HNC 
and ^''CO linewidth distributions are similar, implying that these 
lines are emitted from gas with similar dynamics. In contrast, the 
N2H"'" distribution is signifficantly different, yielding probabilities 
of less than 0.0001 per cent in each case. 

For a temperature of 58 K (see Section l4Tt the median intrin- 
sic thermal linewidths for the molecules are as follows: CH3OH: 
0.29kms~\ "CO: 0.31kms"\ HNC: 0.32kms"^ and N2H+: 
0.31 km s^^. These values represent only 7 to 10 per cent of the 
measured linewidths. Non-thermal line broadening may arise from 
bulk motions, rotation or turbulence, which are difficult to distin- 
guise at the resolutions of these observations. We can conclude that 
the N2H"'" emission emanates from a more quiescent region than 
the other molecules. 



4 DERIVED PARAMETERS 

We present the general physical parameters of the sources and the 
excitation temperatures, column densities and optical depths de- 
rived from the line profiles. The frequencies and electronic con- 
stants used in the calculations are collected in Table [T] 



4.1 Gas mass and virial mass 

Gas masses have been calculated from the 1 .2-mm continuum flux 
density, taken from the work of iHill et ai] l l2005l) . For a 1.2-mm 
integrated continu um flux density Si^ the mass of gas is given by 
lHildebrandl ( [T983h 

Mgas = T, (2) 

±id tiv{^ dust ) 

where D is the distance to the source in metres, is the mass 
absorption coefficient per unit mass of dust, B„ is the Planck func- 
tion for a blackbody of temperature Tdust, and Rd is the dust to 
gas mass ratio. We have adopted a value of 0.1 m^kg^^ for 
at 1.2-mm and have assume d a dust to gas mass ratio of 0.01 
jOssenkopf & Hennind ri994b . The distance has been estimated 
from the Vlsr of th e source using the Galactic rotation curve of 
lBrand&Blit3 ( ll993l) (see Paper 1). Values for Tdust were taken 
from the cold component of the greybody-fit to the spectral en- 
ergy distribution (SED) presented in Paper 1. For sources where we 
had insufficient data to perform a fit we adopted the average value 
of Tdust = 58 K, determined from 65 good fits. We note that these 
temperatures are not well constrained by the sparsely sampled SED 
and are generally higher than the values of 20 - 30 K adopted in the 
literature. For dust temperatures of 20 K and 30 K the gas-mass de- 
rived here will be under-estimated on average by factors of 3.6 and 
2.1, respectively. Values for Mgas were derived for 66 sources and 
are recorded in column six of Table 7. 

We estimated virial masses from the CH3OH, N2H+, 
^•^CO and HN C line-profiles following the approach by 
iMacLaren et al.l (1988). In the simplest case, neglecting support by 
magnetic fields or internal heating sources, the total mass of a sim- 
ple spherical system is given by 
2 

Mvir = ^ = 126rAv2 ... in(MQ), (3) 
where a is the full 3-dimensional velocity dispersion, r is the 
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dust-radius of the cloud in pc, Av is tlie FWHM linewidth of 
the molecular line in kms~^ and G is the gravitational constant 
in Nm^kg"^. The virial mass is the minimum mass required in 
order for a cloud to be gravitationally bound, i.e. the cloud is 
bound if Mgas/Mvir > 1. If magnetic field support is important 
the virial mass may b e overestimated by up to a factor of two 
jMacLaren et alj[l988l) . Observationally, the linewidths may be ar- 
tificially broadened due to blending of multiple components along 
the line of sight, or in optically thick lines, due to radiative transfer 
effects. Enhanced linewidths will cause us to overestimate the virial 
mass by an unknown amount and the values quoted here should 
be considered upper limits. Columns 7-10 of Table 7 present the 
virial masses derived from the four molecules. 

The distributions of Mgas and Mvir are illustrated in Figure|4] 
Optically thin N2H^ and H^^CO"'" produce similar estimates of the 
mean virial mass at 573 M0 and 905 M©, respectively. ^'^CO and 
CH3OH have the broadest linewidths, most likely due to optical 
depth effects or low velocity outflows, leading to mean mass esti- 
mates of 1293 Mq and 1305 M0, respectively. The values derived 
for Mgas and Mvir are consistent with having the same magnitude 
in the majority of clumps. The large uncertainties associated with 
the calculations prevent the interpretation of individual objects as 
gravitationally bound. 



4.2 H2 column and volume density 

We estimated the H2 column density, Nhj , from the gas-mass using 
the following relation: 

M^as 



2MpfHc A' 



(4) 



where Mp — 1.673 x 10~^^kg is the mass of a proton, fno = 
1.36 is a correction factor to a ccount for Helium present in the 
interstellar medium ( lAllenl f 1973) and A is the projected surface 
area of the mm-continuum emission. Assuming a cylindrical source 
with uniform density, the projected area is given by 

A = 7rr^D^ (5) 

where r is the angular radius of the source in radians and D is the 
distance to the source in metres (1 pc = 3.08568 x 10^" m). In a sim- 
ilar manner the volume density can be estimated by replacing 
A in Equation |4] by V the projected volume of emitting gas 



V : 



— TT r D . 

3 



(6) 



Equation |6] assumes emission from a uniformly dense sphere of 
angular radius r. The radius of the dust continuum emission was 
measured directly from th e 1.2-mm continuum images (provided 
courtesy of lHill et alj200A by fitting a Gaussian to the azimuthally 
averaged source profiles. We find values for the full-width half- 
maximum vary between 99 " and 24 " (the limiting resolution of 
the SIMBA bolometer), and have a mean value of 34". Derived 
values of Nhj and are presented in columns four and five of 
Table 7 and have means of 6.1 x 10^^ cm"^ and 4.8x 10'* cm"^, 
respectively. 

We urge caution in interpreting any value of the H2 column 
density derived in this manner as the assumptions introduce large 
unknown errors. Firstly, the mass absorption coefficient and dust 
to gas ratio used to calculate Mgas in Equation [2] may vary con- 
siderably from source to source, leading to a corresponding error 
in the H2 column density. The spatial distribution of II2 is a fur- 
ther unknown and additional errors are introduced by assuming a 



uniform column density over a spherical projected volume. For the 
six sources at the limit of SIMBA resolution (24 ") derived values 
should be considered beam averaged lower limits. Also see Sec- 
tionl6.3lfor a discussion on beam dilution. 



4.3 CH3OH rotational temperature and column density 

CH3OH is a slightly asymmetric rotor, resembling a symmetric top 
except for the 0-H group, which is angled with respect to the prin- 
cipal axis. Internal hindered motion coupled with the rotation of 
the molecule results in a complicated rotational spectrum, organ- 
ised into non-degenerate A and two-fold degenerate E levels, cor- 
responding to different torsional symmetry. The A and E states have 
nuclear statistical weigh ts of 2 and 1 , respectively, thus their overall 
degeneracy is the same teockelee-Morvan et al .Ill994h . 

A rotational temperature Trot and column density N may be 
estimated from the ratios of the thermal CH3OH lines based on the 
assumptions of low optical depths and that the background temper- 
ature (Tb ~ 2.7 K) is small compared to the brightness tempera- 
ture of the line. Under conditions of Local Thermodynamic Equi- 
librium all the observed energy levels can be described by a sin- 
gle excitation temperature Tex, defined by th e Boltzmann relation. 
Using the relation for rotational temperature ^Goldsmith & LangeJ 
\199% . 

N \ 



Q(T) 



kTr, 



(7) 



a rotation diagram may be drawn by plotting ln(Nu/gu) versus 
Eu/k. The rotation temperature and column density in the upper 
level Nu are found from the slope and y-axis intercept, respectively, 
of a straight line fit to the data. The total column density N may then 
be found assuming local thermal equilibriu m and a partition func- 
tion Q (Trot), which for CH3OH is given by jTownes & Schawlowl 
[T95I 



Q(T) = ^(2J + 1) e"'^'/'"'^ « 1.2327 TJ, 



(8) 



In Equation[8]we have appr oximated the partitio n function by a fit 
to the discrete values in the lPickett et al.l (Il9981 spectral line cat- 
alogue. Note that values for Eu are referenced to different ground 
state energy levels in A and E-type CH3OH. The ground state for 
A-type is the Jk ~ Op level, while for E-type CH3OH is the 
Jk = 1-1 level jMenten et al.lll988[) . This is commonly over- 
looked in the literature as all values quoted in' Pickett et alj ( IT998I) 
are referenced to the Jk = level, resulting in incorrect rotation 
diagrams. The correct values of Eu/k are cited in Table[T] 

An example of a CH3OH rotation diagram is shown in Fig- 
ure |5] (/e/r). The data point at Eu/k~7K corresponds to the A+ 
line, while the other three points correspond to E transitions. In 
general, where we detect three E-type lines, they are well fit with a 
straight line, implying that our assumption of optically thin condi- 
tions is valid. 

A and E type methanol are oblate and prolate forms, respec- 
tively, of the assymetric top and only inter-convert on timescales 
of ~ 10,000 years. It is not clear that the A/E ratio can be assumed 
equal to one and a rotation diagram including both A and E type 
transitions may not be valid. In Figure |5] (left) the solid line indi- 
cates a fit to the E-type transitions only, while the dotted line is a fit 
to both the A and E-type. Rotation temperatures derived from only 
the E-transitions have a mean value of 6.67 K, similar to the mean 
temperature derived from all four transitions, which is 6.8 K. Rota- 
tional temperatures and column densities derived using both fitting 
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Figure 4. Histograms showing the distributions of the virial mass es derived from N2H+, h13C0+, "CO and CH3OH, and the gas masses derived from 
the 1.2-mm continuum emission from the work of iHill et"al] i2005l) . Virial masses should be considered upper limits. The plot on the bottom right shows the 
distribution of the Mgas/Mvir ratio derived from N2H+. The values of Mgaa and Mvir are consistent with having the same order of magnitude, however, 

Mvir, N2H+ < Mvir, others aS Av^ jj+ is < AVothcrs- 




Figure S. Qeji) Rotational diagram of CH3OH for the source G10.32— 0.16, made using the 96.7GHz lines. The point at Eu 7K corresponds to the 
2o ~> lo A+ transition. The three other points are due to transitions of E-type methanol. The solid line and values for Trot and N correspond to a fit through 
the E transitions only. The dashed line represents a fit to all transitions, {right) Distributions of rotational temperatures derived from CH3OH using both A and 
E-type transitions (top), and E-type transitions only (bottom). 
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standard formula teohlfs & Wilsonll2004h . 
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Figure 6. Distribution of optical-depths derived from the ratios of the hy- 
perfine groups of N2H+ . The mean is 0.35 ± 0.42. 



methods are presented in Table 8 and the distributions of rotational 
temperatures are illustrated in Figure|5](ng/iO. 

Column densities are averaged over the telescope beam and 
are likely underestimated. If the bulk of the CH3OH emission stems 
from within a hot core then the typical size of the emitting region 
will be -0.1 pc (Ri 20, OOP AU. lKurtz et alj|200(]h . Columns three 
and four of Table 8 present the angular size and beam dilution factor 
calculated using this assumption. Columns seven and ten present 
the column densities corrected for this assumed dilution factor. 

On inspection of the bulk of the rotation diagrams (available as 
an online supplement) we see a pattern emerge. The A-type transi- 
tion has a consistently higher Nu/gu compared to a straight line 
fit through the E-type transitions. This indicates a greater abun- 
dance of A-type CH3OH and an A/E ratio > 1. Considering only 
sources in which all three E-type lines were detected (36 sources), 
and assuming equal excitation temperatures in A and E lines, we 
derive A/E abundance ratios ranging from 0.47 to 1 .80, with a mean 
of 1.30. These values are comparable to v alues previously derived 
in the literature, e.g. jMenten et al.l ( Il988l) who calculated A/E ra- 
tios between 1.3 and 2.0, using the same transitions, towards warm 
clumps in Orion KL. 



4.4 N2H+ optical depth and column density 

Optical depths were derived from the blended hyperfine compo- 
nents of N2H~'", which assume set ratios under optically thin condi- 
tions (see Section l3Tt . Assuming the linewidths of the individual 
hyperfine components are all equal, the integrated intensities of the 
three blended groups should be in the ratio of 1 : 5 : 2 under opti- 
cally thin conditions. The optical depth is calculated from the ratio 
of integrated intensities ( J Tb dv) of any two groups using the fol- 
lowing equation: 



/ TB,idv _ l-e- 



(9) 



where ' a' is the expe c ted ra tio of tiIt\ under optically thin con- 
ditions. ICaselli et alj i 19951) report anomalous excitation of the 
Fi,F = l,0— >1,1 and 1,2 1,2 components (in our group-3), 
so we determine the optical depth solely from the intensity ratio 
of group- l/group-2. We find that 97 percent of sources have optical 
depths below 1 . Figure|6]plots the distribution of the derived optical 
depths, which range from 0.035 to 2.43, with a mean of 0.35 and 
a standard-deviation of 0.42. Values for the optical depth towards 
individual sources are presented in column two of Table 9. 

Column densities of N2H^ can be found from the following 



N 



SkTTiv^ / Tb dv 



^^e^"/"^"^, (10) 



Auihc3 V 1 - e 

where Q(Tcx) is the partition function, e xtrapolated from a fit to the 
values given by the jPickett et al.lll998l) spectral line catalogue: 

Q(T) =4.198Tex. (11) 

Equation[TO]requires an estimate of the excitation temperature Tox- 
Unfortunately, under optically thin conditions, the excitation tem- 
perature cannot be determined from the N2H+ line profiles, so we 
assume a constant value of 10 K. For excitation temperatures of 
5 K and 20 K the total column density should be multiplied by 0.8 
and 1.6, respectively. Beam averaged column densities derived us- 
ing Equation [TO] are presented in column five of Table 9. In the 
literature, high-resolution maps show that N2H^ tr aces the spa- 
tial ex tent of the dust closely (e.g. lCaselli et a l. 2002: Pirogov et akl 
l2003h . We have estimated t he average FWHM of the dust emission 
for all sources observed bv lffill etal.1 ( |2005|) (see section |42] | and 
use these values to correct for beam dilution. Total column densi- 
ties corrected by a factor (6'dust/^bcam)'^ are presented in column 
six of Table 9. 

4.5 "CO, HNC and HCN column density 

The HNC, HCN and ^"^CO profiles frequently display asymmetries 
which may indicate that they are optically thick in general. Lacking 
an estimate for both the excitation temperature and the optical depth 
we are unable to calculate reliable column densities. 

4.6 Abundances 

The relative abundance X between two species may be found di- 
rectly from the ratio of their volume densities. Assuming both 
molecules occupy the same volume of space, this is simply equiva- 
lent to the ratio of their column densities. 

ni _ Ni 
n2 N2 



X : 



(12) 



If, however, one or both species subtend angles smaller than the 
beam, we must correct for their relative volume filling factor and 
Equationll2lbecomes 

3 



X=^ 



N2 V^*! 



(13) 



assuming spherically symmetric emission with characteristic sizes 
Q\ and Q2 ■ For the single-position observations presented here, the 
volume filling factors are complete unknowns and constitute the 
greatest source of errors in the abundance. Additional errors stem 
from calibration of the brightness scale (~ 30 per cent) and the 
pointing accuracy of the telescope (~ 20 per cent). Abundances rel- 
ative to H2 depend also on the assumptions made in calculating the 
gas mass from the 1.2-mm continuimi emission and caimot be de- 
termined reliably. 



5 LUMINOSITIES AND MASSES 

The relationship between mass and luminosity provides clues to 
the star formation efficiency, stellar content and evolutionary state 
of the clumps. Figure |7] {lefi) is an expanded version of Fig. 2 in 
ICesaronil ^20051) . showing the relationship between luminosity and 
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Figure 7. (/e/i") Plot of luminosity versus dust + gas mass for 55 sources whose properties could be determined (after Fig. 2. ICesaro nil2005h . UCHll regions are 
represented by solid triangles and isolated masers by crosses. The thick solid line traces the luminosity expected if Mgas = Mstars in the clumps, and the stellar 
population is concentrated in a minimum number of stars with masses ^100 Mq. The dashed line is the luminosity expected from the most massive star in a 
Mgas cluster with stellar masses distributed according a Salpeter IMF (a = 2.35), and a fundamental upper stellar mass Umit of 150 Mq ( Weidner & Kroug^ 
|2004) . We have assumed a mass-luminosity relation of L oc M"^. The dash-dot line traces the expected luminosity from the same IMF assuming L oc 
and the shaded region indicates intermediate values. The observed mass-luminosity relation is fit with a single power law of index 0.68 ± 0.09 (thin solid line). 
(right) Histograms of the distribution in the Mgas/Lstars values for UCHll regions and isolated masers. A KS-test cannot distinguish between them. 



gas mass for the 55 clumps whose properties could be determined. 
The luminosity has been measured from a two-component grey- 
body fit to the spectral energy distribution and is directly related to 
the stellar content of the clump. Sub-millimetre emission from as- 
sumed optically thin dust ( Hill et al.i,2005i) has been used to derive 
the total gas-mass of the clumps. In the plot, the thick line marks 
the expected luminosity from a minimum number of stars of mass 
M* <100Mq, assuming Mgas = Mstars- Cesaroni (2005) divide 
a sample of 21 known hot core flinto 'h ght' and 'heavy' categories, 
based on their position above or below this division. 'Light' cores 
have gas masses below ~100Mq and fall above the thick line. 
Their high luminosity-to-mass ratio implies that the mass of their 
stellar content is greater than that of the gas. 'Heavy' cores have 
gas masses above ~ 100 Mq and fall below the thick line. The mass 
of their stellar content is comparable to, or less than, the mass of 
the gas. ICesaronil j2005h have interpreted the Lstars/Mgas ratio in 
the 'light' cores as suggesting only a single embedded massive star 
exists in these clumps. Conversely the 'heavy' clumps are likely 
to contain multiple massive stars which have formed in clustered 
mode. We do not find any 'light' cores in our sample. 

An alternative explanation for high Lstars/Mgas ratios in the 
'light' sources might derive from the age of the star forming re- 
gions. A massive young stellar object at a late stage of evolution 
could be expected to have dispersed a large fraction of its molec- 
ular gas, through the action of stellar winds, intense UV-radiation 
and outflows from young OB stars. This would result in a decrease 
in the measured gas-mass without a corresponding decrease in the 
luminosity. In Figure |7] the points would move horizontally into 
the 'light' portion of the plot. We note that Cesaroni's light cores 
include Orion-KL, and IRAS 20126-1-4014, neither of which are be- 
lieved to be at a late stage of evolution. Orion is the closest massive 



■* We have the following sources in common with Cesaroni: G5. 89— 0.39, 
G9.62-l-0.19, G10.47-l-0.03, GlO.62-0.38, G29.96-0.02 and G3 1.41-1-0.31. 



star forming region and is not clear if it is representative of massive 
star-forming regions in general. The mass of IRAS 20126-1-4014 
(7±3M0) has been measured assum ing Keplerian rotation of a 
molecular disk jCesaroni et Zl bOO^) and relates to the inner ob- 
ject only. On larger scales the gas- mass has been measu red to be 
several hundred solar masses (e.g. lEstalellaetalJll993l) . putting 
IRAS 20126-1-4014 in the 'heavy' category. We find that all of our 
maser selected sample have gas masses above ~30Mq and all 
fall into the 'heavy' category. If the division between 'light' and 
'heavy' cores is an indication of age, then all of our sample are at a 
relatively early stage of evolution, i.e., the powering massive young 
stellar objects have not yet dispersed their natal molecular clouds. 

In Figure |7] the shaded region between the dashed lines rep- 
resents the range of luminosities expected from the most massive 
star in a cluster where Mstars = Mgas, if the initial mass func- 
tion (IMF) follows a Salpeter power law (q — 2.35) and as- 
suming a fundamental upper stellar mass limit of M, = 150 Mq 
(Weidner & Kroupa 2004). The luminosity of the cluster will be 
dominated by the most massive star so this is a good approxima- 
tion to the expected total luminosity. The range in expected lumi- 
nosities derives from the uncertainty in the mass-luminosity rela- 
tion for massive young stellar objects. This may be represented 
by a power law of the form L* oc M'l, where 7 ranges between 
3 and 4 in the literature for OB-stars (e.g.. lLadalll999l : IChurchwel]| 
l2002i) . Our maser selected sample exhibit a rough correlation be- 
tween mass and luminosity, which may be fit with a single power 
law with index 7 = 0.68 ± 0.08. This is a much shallower rise 
than observed in evolved stars, however, the luminosities of the 
most massive clusters are consistent with the 150 Mq stellar cutoff 
jWeidner & KrouDj|200i) . even in the more conservative limit of 
oc Ml 

The derived gas-mass is highly dependent on the temperature 
in Equation[2l The shallow power law might be explained if the av- 
erage dust temperature increases with luminosity and this change is 
not reflected in the sparsely sampled greybody fits to the SED, from 
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Figure 8. (left) Histograms showing the distributions of dust temperature taken from the cold component of the greybody fit to the SED. The top panel shows 
UCHll regions, which have a mean Tcoid "f 62.4 K. while the bottom panel shows isolated masers, which have a mean T^oid of 58.1 K. (right) Histograms 
showing the distributions of CH3OH rotational temperature for UCHll regions (top) and isolated masers (bottom). The mean values are 9.9 K for UCHll regions 
and 6.1 K for isolated masers. A KS-test establishes that the CH3OH rotational temperature distributions are significantly different, returning a probability of 
0.002 per cent. 



which we take our temperature estimates. If we assume the mass es- 
timates are correct to within a factor of a few, then the main resuh is 
that the lower mass clumps are over-luminous and the higher mass 
clumps are under-luminous. 

5.1 UCHii regions: Luminosity or Age? 

UCHii regions (triangles) are in general more luminous and more 
massive than the isolated masers (crosses) in the sample. Individ- 
ual power law fits to the isolated masers and UCHii sub-samples 
yield mass-luminosity relations of L^oi | masers oc ]V[0-59±o.i3 
Lboi|uCHii ]y[O C4±o.i8^ respectively. These fits are consistent 
with the fit to the combined sample, within the errors, and a single 
power law can describe both distributions. Figure[7](ng/!0 is a his- 
togram of the Mgas/Lstars ratio for UCHii regions (top) and isolated 
masers (bottom). As expected from the mass-luminosity plot, the 
distributions are similar and a KS-test cannot distinguish between 
them, yielding a 76 per cent probability of the two sub-samples be- 
ing drawn from the same population. This result implies that for the 
same masses, the two samples have the same luminosities. Early in 
the evolution of a massive young stellar object we would expect its 
luminosity to increase as it accretes matter and more mass is trans- 
fered from the envelope to the star. At the same time the overall 
mass of gas and dust in the clump is expected to remain relatively 
constant, or indeed decrease, as winds from the young OB-cluster 
break up the cloud. By this rationale, if our sample of UCHii re- 
gions are older and more evolved, we would expect them to have 
a lower Mgas/Lstars ratio than the isolated masers. There are sev- 
eral possible reasons why we do not see a difference between the 
mass-luminosity distributions. 

It is possible that systematic errors (which we discuss be- 
low) associated with our assumptions cause us to overestimate the 
mass of the clouds hosting the UCHii regions relative to the clouds 
hosting the isolated masers. For there to be no difference between 
the mean masses of the sub-samples, the isolated masers would 
need to have a dust temperature on average 36 K lower than the 
UCHii regions. Figure [8] (/e/r) illustrates the distributions of Tcoid 
from the SED fits, for UCHii regions (top) and isolated masers (bot- 
tom). While the UCHii regions clearly have temperatures which are 



skewed to the high end, the difference in the means is only 4.3 K. 
The rotational temperatures calculated from CH3CN are likewise 
similar, however, we cannot rule out different average dust tem- 
peratures on the basis that CH3CN emission likely probes only 
the embedded 'hot core' and does not reflect the conditions in the 
extended clump. Thermal CH3OH emission is thought to derive 
from a more extended envelope and rotational temperatures calcu- 
lated from this molecule appear to support very different tempera- 
tures in the two groups. Figure [8] (n^A/) shows histograms of the 
CH3OH rotational temperature for the UCHii regions (top) and iso- 
lated masers (bottom). A KS-test returns a probability of 0.002 per 
cent of the distributions deriving from the same parent. We note 
that the mean values of 9.9 K for the UCHii regions and 6. 1 K for 
the isolated masers do not reflect the kineti c temperature as the 
CH3O H is likely sub-thermally excited (see [Goldsmith & LangeJ 
1 19991) . 

As an alternate estimate of the mass we derived a 'LTE-mass' 
from the column density of N2H+. In our calculations we as- 
sume (i) the N2H^ emission closely traces the 1.2-mm dust con- 
tinuum and has a flat density distribution within the FWHM of the 
clum p, (ii) a fixed [H 2]/[N2H+] abundance ratio of 1.7 x 10^"' 
(Piro govetal]|2003l) . The mass of gas may then be estimated from 
the column density using Equation |4] Figure [9] (/e/f) is is a mass- 
luminosity plot similar to Figure|7] except that the mass plotted on 
the X-axis is the LTE-mass calculated from Although the 

scatter is greater, we clearly see the masses of the UCHii regions 
and isolated masers are still weakly offset, with the UCHii regions 
having greater masses on average, as in Figure[7l The distributions 
of the Mgas/Lstars ratios are illustrated in Figure [9] (ng/zf). A KS- 
test returns a probability of 50 per cent that both distributions are 
derived from the same parent, confirming the original result. 

Another explanation for the offset in the masses derives from 
the sensitivity limit on the radio-survey. In the preceeding paper we 
assumed the UCHii regions were older and more evolved than the 
isolated masers. We argued that the original radio continuum sur- 
vey was sensitive enough to detect UCHii regions associated with 
the less luminous isolated masers (see Paper 1). However, if a sen- 
sativity bias existed it would manifest as a mass and luminosity off- 
set between the radio-loud and radio-quiet, i.e, the UCHii regions 
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Figure 9. As Figure [7] except in this case the mass is the LTE-mass derived from the column density of N2H+ assuming a fixed abundance ratio of 
1.7 X 10""' and that the N2H+ emission closely follows the spacial distribution of the 1.2-mm continuum emission. 



would be more massive and more luminous than isolated masers. In 
other words, some of the lower Lstars isolated masers may indeed 
be associated with UCHii regions which fall below our detection 
limit. In Paper I we showed that most sources are luminous enough 
to have detectable UCHii regions, assuming up t o 90 per cent of th e 
Lyman continuum photons are absorbed by dust jKurtzetal.ll994h . 
If the disparity in luminosity is due to a real sensitivity limit, this 
would imply that, on average, at least 99 per cent of the ionising 
photons are absorbed by dust before contributing to the creation of 
a UCHii region. We believe this is unlikely in the majority of cases. 

Some of the isolated maser sites may also be associated with 
hyper-compact Hit (HCHii) regions. This recently dis covered class 
of H ii region is typically less than O.OI pc across jAvalos et al.l 
I2OO6 ) and exhibits a free-free emission spectrum with a turnover 
at higher frequencies than UCHii regions, typically 20 GHz. HCHii 
regions are expected to be optically thick at 8.6 GHz due to the 
squared dependac e of the flux density wi th frequency and will be 
difficult to detect. iLongmore et al.l ( l2007h had used the Australia 
Telescope Compact Array to search for 22 GHz emission towards 
24 of the objects presented here. Nine of these sources fall into the 
radio-quiet category and of these seven exhibit compact (30 -^3") 
continuum emission. This result strongly suggests that some of our 
radio quiet sources are associated with HCHii regions, or indeed 
UCHii regions undetected at 8.6 GHz. However, because HCHii re- 
gions are thought to be the precursors of UCHii regions, we are still 
be selecting for two different evolutionary states. 

The most obvious explanation for the offset is due to the 
timescales on which we expect massive young stellar objects to 
progress to develop a detectable UCHii region. In essence, clumps 
containing very high mass objects will evolve very rapidly and will 
spend only a short time in hot core and earlier phases. Because 
of the short timescales involved, very massive hot cores will also 
likely be deeply embedded. The first indication of star-formation 
activity may be when the expanding Hit region begins to dis- 
perse the surrounding molecular gas. Anecdotal evidence support- 
ing this explanation is already emerging from uniform searches for 
MYSOs, such as the Red MSX Source (RMS) survey jHoare et al] 
l2004lMottranill2008h . The less evolved MYSOs identified in the 
survey have a significantly smaller upper-mass-cutoff compared to 



the UCHii regions detected. No MYSO counterparts are found to 
the most massive UCHii regions, implying that the UCHii phase is 
the first indication of star-formation in the most massive clumps. 

We believe differences between the radio-bright and radio- 
quiet samples may be attributed to evolution, with the caveat that 
the highest luminosity objects do not spend significant time in an 
'isolated maser ' phase. 



6 DETAILED COMPARISON BETWEEN SUB-SAMPLES 

The aim of this molecular survey is to investigate the utility of 
molecular lines in the 3-mm band to act as 'molecular clocks'. 
In this section we consider the differences between sub-samples 
divided from the source list based on associations with tracers of 
potentially different evolutionary states. 

In Paper 1 we examined the link between hot cores, 6.67 GHz 
CH3OH masers and UCHii regions. We found that CH3CN was 
brighter and more commonly detected towards masers associated 
with UCHii regions, than towards 'isolated' masers. Assymetries 
in the HCO^ versus H^"^CO^ line profiles were interpreted as in- 
ward (blue skewed) or outward (red skewed) gas motions. We also 
examined the associations with mid-infrared emission in the Mid- 
course Space Experiment (MSX) images and found 1 1 dark clouds 
(seen in absorption against the 8 p.m Galactic background) and 67 
in emission. Many of the sources also exhibited high velocity line 
wings, indicative of outflows. Based on these results we separate 
the source list into sub-samples based on the following criteria: 

• Presence or absence of radio emission from a UCHii region. 

• CH3CN detected versus no CH3CN detected. 

• MSX-dark (seen in absorption at 8 fim) versus MSX-bright. 

• Presence or absence of high velocity line- wings in HCO^, 
HNC or ^^CO. 

• Presence or absence of blue-skewed HCO"'" line profiles, indi- 
cating infall. 

In the following analysis we contrast the sub-samples by looking 
for differences in the median properties. A KS-test is used to assess 
the similarity of the two distributions and trends are found based 
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Figure 10. The bars in this horizontal histogram represent the difference 
between the median values of two distributions, normalised to the standard 
deviation of the combined sample. Here, we compare spectral line luminos- 
ity in sources with and without UCHll regions. A positive value means the 
line is more luminous towards UCHll regions. To assess the significance of 
the differences, we compute the KS-statistic (left column) and associated 
percentage probability (right column). As a visual aid, the greyscale of the 
bars is scaled to the probability. Probabilities smaller than approximately 
two per cent (dark bars) indicate significant differences. 
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Figure 11. Differences in the median Hnewidth, for sources with and with- 
out an associated UCHll region. The linewidths were measured in two ways: 
by fitting a Gaussian to the line profile and measuring the full-width at half- 
maximum, or by integrating under the line and computing the equivalent 
width (HCN only). Other details are as in Figure fTol 
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Figure 13. Differences in the median line-line column density for sources 
with and without an associated UCHll region. 
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Figure 12. Differences in the median line-line intensity ratio for sources 
with and without an associated UCHll region. For molecules with multiple 
transitions (e.g., CH3CN, CH3OH) we have summed the intensity of all 
measured lines. 
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on large differences in the median values, coupled with a low prob- 
ability that the sub-populations are drawn from the same parent 
distribution. 

Initially, we compare the spectral line luminosities calculated 
from the integrated line intensity multiplied by the square of the 
distance to the source. If the distance is well known then com- 
parisons of line luminosities will potentially reveal intrinsic dif- 
ferences between the sub-samples. However, the error in the dis- 
tance for this sample is relatively large (±0.65 kpc ~ 20 per cent), 
and increases drastically towards Galactic longitudes < 5°, where 
the model Galactic rotation curve begins to break down. Only the 
strongest trends will be seen using such a comparison. 

We also compare the distance-independent linewidths. 
Linewidths in massive star forming regions are generally much 
greater than the thermal width and are indicators of the level of 
turbulence present in the gas. Turbulence arises from bulk mo- 
tions of gas, is often associated with energetic phenomena, such 
as outflows, and is correlated with star formation activity (e.g. 
iDobbs et al.ll2005h . 

The ratio of two spectral line intensities is unaffected by the 
distance and constitutes a much better probe of the physical and 
chemical conditions. The dominant errors are due to the relative 
accuracy with which the data has been calibrated and the unknown 
beam-filling factors. The results presented here are based on data 
taken over four years and the calibration error in the data has been 
determined to be better than 30 per cent in most cases. The uncer- 
tainty due to the beam filling factors is given by the ratio of the solid 
angles over which the two molecules emit, {Q.\/Q.2), and may be 
large for dense-gas tracers which subtend angles smaller than the 
beam. The trends underlying the data should still be visible given 
the sample size. 

Ideally, in a search for molecular clocks the abundances 
should be compared directly. Unfortunately, the small number of 
transitions observed in this work make it difficult to constrain the 
column densities accurately. Two or more transitions of a single 
molecule are needed to solve for the optical depth, excitation tem- 
perature and column density. We have been able to determine the 
column densities of thermal CH3CN and CH3OH through the ro- 
tational diagram method, and by assuming an excitation tempera- 
ture or [^^C]/[^''C] abundance ratio, we have calculated the col- 
umn densities of N2H+ and HCO^. Here we compare the beam- 
averaged column densities of these molecules directly to find the 
abundance ratios. 

Lastly, we contrast the differences in a suite of measured 
physical parameters such as bolometric luminosity, gas-mass, 
rotational-temperature and mid-infrared flux density, drawn from 
this and other work. 

The important differences are sunmiarised in section |6!4l 



6.1 Radio-loud versus radio-quiet sources 

Figure [TO] illustrates the difference in line luminosity between 
sources with and without 8 GHz radio continuum emission (from 
UCHii regions) in the Mopra beam: the radio-loud and radio-quiet 
sub-samples, respectively. As with CH3CN and HCO^in Paper 1, 
we find that the luminosity of the remaining transitions is enhanced 
towards the radio-loud sub-sample. The median line luminosity is 
greater by 0.65-(7, on average, where o is the standard deviation 
of the combined sample. Thermal CH3OH stands out as having the 
most significantly different distribution and the largest difference in 
line-luminosity, followed by H^3C0+, CH3CN, HNC and N2H+. 
We find the median linewidths of all species are broader to- 
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Figure 14. Differences in other computed and measured parameters, for 
sources with and without an associated UCHll region. 



wards the radio-loud sub-sample (see Figure [TTt. CH3OH exhibits 
the greatest differences, closely followed by N2H"'"and H^^CO^, 
all of which have statistically different distributions according to 
the KS-test. HCO+, HCN, HNC and CH3CN have broader lines 
on average, however the distributions are not significantly different. 
The differences in the median linewidth between the radio-loud and 
-quiet subsample range from 0.37 kms~^ to 1.62kms~^, greater 
than can be attributed to thermal broadening and most likely due 
to the greater turbulent energy in the gas associated with UCHii re- 
gions. 

The differences found between the linewidths and line- 
luminosities are echoed in a comparison of the raw integrated 
line intensities. Figure [T2| shows the difference in the median val- 
ues of the 28 intensity ratios between eight molecules. The me- 
dian intensity of CH3OH and, to a lesser extent, CH3CN, are 
clearly enhanced with respect to all other molecules. Statistically, 
the most significant difference in intensity ratio occurs between 
CH3OH and other molecules, specifically N2H+ and ^^CO. The 
[CH3CN/^^CO] and [H^^CO+Z^^CO] ratios also stand out as be- 
ing significantly greater towards the radio-loud sub-sample. 

Figure [T3] illustrates a comparison of the median beam- 
averaged column densities of four molecules: CH3CN, CH3OH, 
HCO^ and N2H"''. Again, thermal CH3OH stands out as having 
greater a column density, and hence a higher abundance, towards 
the radio-loud sub-sample. 



Physical and chemical conditions in methanol maser selected hot-cores and UCHII regions 13 



DIFFERENCE IN LINE LUMINOSITY 

[HNC] X g 
[HCO*] X 

[HCN] X 
[h"C0+] X 
[CH3OH] X 

["cO] X 

No CH3CN , 



KS 


%-PROB 


0.23 


39.4 


0.26 


23.8 


0.34 


5.1 


0.32 


7.8 


0.19 


62.7 


0.24 


42.3 


0.26 


23.8 


CH3CN , 





-0.1 0.1 

(Median. —Medlan2) / Sigma 



N^H* FWHM 
HNC FWHM 
HCO* FWHM 
HCN E-Width 
h"cO+ FWHM 
CH,OH FWHM 



DIFFERENCE IN LINEWIDTH 

P 



KS %-PROB 



0.18 67.3 



0.17 74.1 





■ 


1 0.38 


2.5 






0.14 


91 .7 




□ 


0.13 


95.8 




1 


1 0.32 


12.8 


■ 




0.31 


1 1.2 


3CN , 




CH3CN , 
... 





-0.5 0.5 

(Median^— Iv1edian2) / Sigma^ 



Figure 15. Differences in tlie median line luminosity, for sources with and 
without CH3CN. 



Figure 16. Differences in the median linewidth, for sources with and with- 
out CH3CN. 
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Figure 17. Differences in the median line-line intensity ratio for sources 
with and without CH3CN. For line with multiple transitions we summed the 
integrated intensity of all measured lines. 
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Figure 18. Differences in the median column densities, for sources with 
and without CH3CN. 



The median differences between other measured and calcu- 
lated properties are shown in Figure [14] The median bolometric 
luminosity and median mass of the radio-loud sub-sample is signif- 
icantly greater than for the isolated masers (see section[5)- The rota- 
tional temperature derived from CH3CN is similar in both groups, 



however, the CH3OH rotational temperature is significantly greater 
in the radio-bright group (also see Figure[8). 
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6.2 Sources with and without CH3CN 

CH3CN was detected towards 58 sources (70 per cent), however the 
detection rate towards the radio-loud sub-sample (19 sources) was 
90 per cent. The enhanced luminosities and masses of our UCHii re- 
gions may bias the results of a comparison between CHsCN-bright 
and CHaCN-dark groups. We therefore exclude any regions with 
radio emission from the comparison. In practise, we find that the 
major results are robust to the presence of UCHii regions. 

The differences in median line luminosity between the two 
groups presented in Figure [TS] are not as pronounced as in the pre- 
vious section, however, all lines are marginally (O.OV-tj) more lu- 
minous towards the CHaCN-bright sub-sample. HCO^ and HCN 
have the most significantly different distributions. We note that the 
HCO'*' and HCN profiles are affected by self absorption, with the 
CHsCN-bright group containing the most extreme examples. The 
significant difference in the HCO"'' luminosity is also not reflected 
in the optically thin H^'^CO^ isotopomer, suggesting that an abun- 
dance difference is unlikely. The greatest difference in the median 
values occurs in "CO, followed by CH3OH and H^^CO+. 

Linewidths are also biased towards being greater in the sample 
with detected CH3CN (see Figure[T6t. HCO^ exhibits the most sig- 
nificant difference, followed by CH3OH, although a KS-test cannot 
distinguish between the populations with confidence. 

When we examine the line intensity ratios in Figure[T7]a sim- 
ilar picture emerges to the case of the UCHii regions. Thermal 
CH3OH stands out, as its intensity relative to all other lines is en- 
hanced in the CHaCN-bright sub-sample. We also see evidence 
for enhanced N2H+ emission. In contrast, ^"'CO emission is rel- 
atively weak towards the CHsCN-bright sub-sample. These ratios 
are highly suggestive of different environmental conditions in the 
dense gas, (traced by N2H^ and CH3OH) compared to the more 
extended gas (traced by ^^CO), and correlate well with the line lu- 
minosities. 

Figure [TS] shows the differences in the median column densi- 
ties of CH3OH, HCO+ and N2H+. Again, the thermal CH3OH col- 
umn is seen to be significantly enhanced towards the CH3CN bright 
group, echoing the result from comparing the raw line ratios. 

The difference in the median values of other derived parame- 
ters are displayed in Figure|20l Sources with detected CH3CN have 
significantly brighter thermal continuum emission at I.2-mm wave- 
lengths. This translates into a mean mass higher by ~300Mq, 
omitting UCHii regions or ~600Mq, including UCHii regions. 
Histograms of the mass distributions for the two groups are plot- 
ted in Figure [T9](/e/0, and show significant differences. All other 
parameters presented in Figure|20lare similar for both sub-samples. 

The division of the sample into populations with and without 
detected CH3CN may be a somewhat artificial distinction. While 
CH3CN emission may be present in more distant sources, it may 
fall below our 2-a detection threshold of 80 mK. Emission over a 
small solid angle will suffer from high beam-dilution factors, fur- 
ther compounding the problem. The distributions of distance in the 
two samples, plotted in Figure[T9] (ng/;f), are similar, except for a 
tail of more distant objects without CH3CN. This difference sug- 
gests that a small fraction of non-detections may be attributable to 
sensitivity limits. 

We also split the sample via other criteria (see the introduction 
to this Section), however, the chemical and physical differences are 
not as significant. The results of these investigations are presented 
as an online supplement. 
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Figure 20. Differences in other computed and measured parameters, for 
sources with and without CH3CN. 



6.3 The effect of beam dilution 

As mentioned in Section l431 the beam filling factors constitute an 
unknown source of error when comparing the intensities or column 
densities of two molecules. If both emitting regions subtend angles 
(9) smaller than the beam, then we must correct for the relative 
volume filling factor, given by (61/62)^- Assuming the emitting re- 
gions have similar characteristic sizes in all sources, the relative 
beam filling factors will average to ~ 1 over the entire sample. 
However, if one species generally fills the beam while the other 
species is generally beam-diluted (i.e., one has constant angular 
size while the other doesn't, then we would expect a 1/distance^ 
dependence in the line intensity ratios). Figure [2T] presents two 
plots of line intensity ratio versus distance. The left panel plots 
the CH3CN/CH3OH ratio versus distance. The emission from both 
molecules is expected to derive from regions smaller than the 
36" Mopra beam, even for nearby sources. We see the line-ratios 
are not correlated with distance, consistent with a constant ratio 
between the angular sizes of the species. The right panel plots 
the CH3CN/N2H^ ratio versus distance. Ass uming N2H~'' emis- 
sion closely follows the dust morphology (e.g.. lCaselli et alj|2002l : 
|Pirogovetal1l2003l) , we have shown in Table 9 that the beam- 
dilution factor is generally close to I. In Figure |2T] (right) the 
dashed line plots the expected effect of CH3CN beam dilution as a 
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Figure 19. (left) Histograms showing tlie distributions of gas masses derived from 1.2-mm continuum flux density for sources with (top) and without (bottom) 
detected CH3CN. Sources in which CH3CN has been detected tend to be more massive, (right) The distribution of distances for sources with (top) and without 
(bottom) detected CH3CN. A tail of more distant objects exists in the CH3CN-dark sample, suggesting that a small fraction of non-detections may be due to 
sensitivity limits. 



source is moved further away. We see the scatter in measured hne- 
ratio increases with distance making it difficult to interpret the data 
and, in addition to beam dilution, chemical and physical differences 
in the sources affect the line intensity ratios. Omitting the UCHii re- 
gions (squares) and considering only isolated masers (crosses), we 
see a possible trend towards decreasing CH3CN/N2H+ ratio with 
increasing distance. This trend does not approach the 1/distance^ 
dependence expected, but it may indicate that some of our compar- 
isons will include beam dilution effects for CH3CN. 

6.4 Interpretation and discussion 

It is clear we are seeing significant differences between sources as- 
sociated with different tracers. In this section we summarise the 
differences between the sub-samples and suggest some possible in- 
terpretations of the results. 

U CHii regions vs isolated masers 

The clearest distinctions exist between isolated CH3OH masers and 
those with UCHii regions. The result found in Paper 1, where the 
CH3CN and HCO^ lines were shown to be more luminous and 
have greater linewidths in the presence of a UCHii region, has been 
extended to the other observed species. When we examine the raw 
line intensity ratios, thermal CH3OH distinguishes itself as being 
relatively intense towards the radio loud sources. This is reflected 
in both the higher median column density and CH3OH rotational 
temperature. We interpret this result as evidence for an enhanced 
abundance of CH3OH in relatively warm gas nearby the UCHii re- 
gions. Significant differences are also found between other tracers 
of high and low density gas. The [CH3CN/"C0], [N2H+/^^CO] 
and [H^^CO+Z^^CO] ratios are enhanced in the radio-loud sub- 
sample, possibly indicating different excitation conditions or abun- 
dances in the 'core' compared to the 'envelope' gas. Interestingly, 
when compared to all other lines, the relative brightness of the hot- 
core tracer CH3CN is not significantly greater, yet we detect it to- 
wards 95 per cent of UCHu regions and only 63 per cent of the 
isolated masers. 

A potential explanation for the median offset in line- 
luminosity and linewidths between the radio -quiet and -loud sam- 
ples may be found in the significantly greater bolometric luminosity 



and gas-mass of the radio-loud sample. Simply put, more luminous 
clumps will likely have a greater flux of IR-photons and higher gas 
temperatures. Under these conditions one would expect to see more 
luminous spectral lines. Linewidth is commonly used as a proxy for 
star formation activity, and hence evolutionary state, as it is largely 
dependent on turbulence in massive star forming regions. However, 
on the scale of molecular clouds (3-20 pc) the lin ewidth has also 
been empirically correlated with the cloud mass J Solomon et al.l 
Il987h . possibly due to multiple (sub)-clumps with different veloc- 
ities within the beam. It is unclear which effect dominates in our 
sample. Figure [22] plots the linewidth N2H+ versus the gas-mass. 

is chosen as it is optically thin and has been shown to trace 
the dense regions of the molecular clouds (e.g., ,Caselli et al.,,2002l) . 
A weak correlation is observed which may account for some of the 
linewidth differences between the groups, however we cannot rule 
out an evolutionary effect. 

That the relative intensity and abundance of CH3OH is greater 
towards the radio-loud sources likely reflects real chemical differ- 
ences between the sub-samples. This may be interpreted as evi- 
dence for the more advanced evolutionary state of the radio-loud 
sample. The CH3OH abundance is predicted to be enhanced by a 
number of mechanisms as a massive young stellar object evolves. 
In the hot-core phase it is known to b e evaporated from th e dust 
grains as the heating source evolves jChamlevet .^119951) . The 
abundance of CH3OH is also predicted to be enhanced in the walls 
of outflow caviti es, where it is spu ttered from the grains by low- 
velocity shocks l lHogerheiidell2005b . 

CHzCN-bright vs CH^CN-dark 

In this and the remaining comparisons we removed the radio-loud 
sources from the list, as their higher intrinsic luminosities may bias 
the results. We find sources with detected CH3CN are also biased 
towards having greater line-luminosities and linewidths, although 
these differences are not as pronounced as with the radio-loud ver- 
sus radio-quiet division. The most striking result is the significant 
tendency for thermal CH3OH to be brighter and have a greater 
column-density in the CH3CN bright group. This is not unexpected 
as the CH3CN and C H3OH abundances ar e predicted to increase 
with time in hot cores jChamlevetaill 19951) . We also see evidence 
for more turbulent and energetic conditions in the dense gas (traced 
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Figure 21. Plots of line-line intensity ratio versus distance, (left) [CH3CN/CH3OH] ratio vs distance, (right) [CH3CN/N2H+] ratio vs distance. Crosses 
mark isolated masers, while squares mark UCHll regions. The dotted line show how the [CH3CN/N2H+] ratio is expected to decrease with distance if the 
CH3CN emission has a fixed physical size and the N2H+ emission always fills the beam. 



by N2H and CH3OH) compared to more extended and diffuse 
gas (traced by ^''CO), implying a more advanced evolutionary state 
in the CH3CN-bright sub-sample. Interestingly, the CH3CN-bright 
sub-sample also has a significantly greater median mass, which 
may bias this result. CH3CN is also not detected in sources with 
distances greater than 10 kpc, suggesting that a small number of 
sources may contain CH3CN below the sensitivity limit of our sur- 
vey. 

MSX-bright vs MSX-dark 

All species in the MSX-dark sample have marginally greater line- 
luminosities and linewidths compared to their MSX-bright coun- 
terparts, with the ^'^CO linewidth being especially notable. This 
is unexpected, as infrared dark clouds have been postulated to be 
the young and q uiescent precursors to hot cores and UCHii regions 
dRathbome etalj.2006) . HCO+ and CH3OH tend to be brighter and 
CH3OH marginally more abundant in the MSX-dark sub-sample, 
although this is not judged to be significant. Other parameters are 
similar in both sub-samples. Given the small number of sources in 
the MSX-bright subsample (9-11 sources), none of the observed 
differences are large enough to draw any definite conclusions from 
the data. 

Line-wings vs no line-wings 

There are no significant differences between the median line- 
luminosities when the sources are divided into sub-samples with 
and without high-velocity line-wings. Differences between the 
HCO"'' linewidth distributions can be attributed to unidentified low- 
velocity linewings. ^"^CO is significantly brighter (relative to other 
lines) in the sources without linewings and CH3CN is biased to- 
wards being brighter and more abundant in the sub-sample with 
linewings. No distinguishing differences are seen in the physical 
parameters or IR-colour ratios. 

Infall profiles vs no infall profiles 

The properties of sources with and without HCO"'' 'infall-profiles' 
are close to identical. Little difference is seen between the median 
line-luminosities, linewidths or line intensity ratios. Biases in the 
H"C0+ and HCO+ line ratios can be attributed to selection ef- 
fects. 

The online material contains additional discusion and figures 
pertaining to the previous three paragraphs. 
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Figure 22. Plot of N2H+ FWHM linewidth versus gas-mass. Squares 
represent UCHll regions and crosses represent isolated masers. 



6.5 Further discussion 

Is the mass of the hot core related to the dust mass? Figure[23lshows 
a plot of CH3CN beam averaged column density versus mass of gas 
derived from 1.2-mm continuum emission. A clear positive correla- 
tion is seen between the two. If we assume the column of CH3CN is 
a proxy for the mass of the hot core, then it follows that the mass 
of the hot core increases with the mass of the extended clump. The 
relationship is well fit by MchsCn (x M^uat and assuming the 
density is constant within the gas, this imp lies that Mhot co re oc 
Mdusf This is consistent with the results of ICesaronil feOOS *! who 
found that the hot core is an continuation of the clump, and there is 
no break in the approximately 1/r^ density power law. Rather, the 
hot core is differentiated by temperature and excitation conditions 
influenced by the central heating source. 



7 SUMMARY AND CONCLUSIONS 

In continuation of the 'Hot Molecular Cores' survey we used the 
Mopra telescope to search for 3-mm transitions of CH3OH (2-1), 
N2H+ (1-0), "CO(l -0), HCN (1 -0) and HNC (1,-0). This is in 
addition to the previously reported measurements of HCO^ (1-0), 
H"C0+ (1-0), and CH3CN(5-4) & (6-5), in Paper 1. Molec- 
ular emission was detected in all but one source (GIO.IO— 0.72), 
with N2H+, ^'^CO, HNC and HCN detected towards all of the re- 
maining 82 targets. Thermal CH3OH emission was detected in 78 
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Figure 23. Plot of CH3CN column density versus gas-mass derived from 
1.2-mm continuum flux density. 



sources (94 per cent). The following conclusions have been drawn 
from the data: 

(i) Virial masses estimated from N2H+ are comparable to or 
lower than the gas masses of the sample. Values for Mgas have 
been calculated from the 1.2-mm dust emission and range from 
30 M0 to 5 X 10^ Mq, while values for Myiriai range from ~40 
to 2.5 X 10^ M© . We note that N2H+ exhibits particularly narrow 
linewidths with an average of 3 .0 km s ~ ^ . 

(ii) H2 volume densities have also been estimated from the sur- 
face brightness and extent of 1.2-mm emission. The mean value is 
5 X 10'' cm~^, however, we note that the assumptions inherent in 
the calculation introduce large unknown uncertainties in values for 
individual sources. 

(iii) Rotational temperatures derived from CH3OH range from 
3.0 K to 14.0 K, with a mean value of 6.67 K. These anoma- 
lously low temperatures strongly suggest that the CH3OH is sub- 
thermally excited. Assvuning a common excitation temperature, we 
find the abundance ratio of A to E-type CH3OH ranges from 0.5 
to 2.8, with a mean of 1.5, consistent with previously published 
values. 

(iv) N2H"'" is found to be optically thin towards most of the sam- 
ple. Likewise the CH3OH rotational diagrams are consistent with 
low optical depths in the J = 2 - 1 transitions. CO, HCN and HNC 
all exhibit line profile asymmetries and may be significantly opti- 
cally thick. 

(v) The luminosity/gas-mass relationship for the sample is 
found to follow a rough power law of the form LocM" ®*, mean- 
ing that the lower mass clumps are over-luminous and the higher 
mass clumps are under-luminous. We find that the gas-mass of the 
clumps is comparable to or greater than the mass of the stellar con- 
tent (implied by the luminosity). We interpret this as an indicator 
of young ages, as the newborn massive stars still retain their natal 
molecular clouds. 

(vi) The radio loud sub-sample (UCHii regions) is on average 
more massive and more luminous than the radio-quiet sub-sample 
(isolated masers). We argue that the the most likely explanation 
for this disparity is due to the comparatively rapid timescale over 
which the most massive objects evolve into Hii regions. The first 
significant tracer of star-formation in these objects may be radio- 
continuum emission from a HCHii or UCHn region. It is also possi- 
ble that some of the isolated maser sample contains hyper-compact 
Hii regions wliich are too optically thick to detect at 8 GHz. In ei- 
ther case there is a clear age difference between the subsamples. 

(vii) We separated the source list into sub-samples associated 



with different tracers of evolution and compared their median phys- 
ical and chemical properties. Our findings are as follows: 

• All spectral lines are brighter and more luminous towards 
the radio loud sample, extending the result found in Paper 1. 

• The thermal CH3OH brightness and abundance is found 
to be enhanced (compared to other species) in radio-loud and 
CHsCN-bright sources. We suggest that in single-dish observa- 
tions thermal CH3OH may be used as a crude molecular clock 
to indicate the evolutionary state of a clump. 

• Given the small sample size, conclusions drawn from a 
comparison of MSX-bright and MSX-dark sources are not sig- 
nificant. It is worth noting, however, that linewidths and line- 
luminosities are generally greater in the infrared dark clouds, in- 
dicating perhaps that when selected using CH3OH masers, dark- 
clouds conceal objects at later stages of evolution. 

• We divided the source list into sub-samples with and with- 
out high velocity linewings and HCO"'"-infall-profiles, but found 
no significant differences in chemical or physical properties. 

(viii) We find that the beam-averaged CH3CN column density is 
correlated with the gas-mass derived from 1.2-mm emission. This 
suggests that the mass of the hot-core is dependent on the mass of 
the clump in which it is embedded. 
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Table 1. Details of observed transitions. 



Species 


Transition 


Frequency 
(GHz) 


(K) 


(s-i) 




Su 


CH3OH " 


2(-i,2) ^ 


96.739390 


4.642 


2.495x10" 


-6 


5 




•^(0,2) (Oil) 


96.741420 


6.963 


3.327x10" 


■6 


5 




2(0,2) ^ 1(0,1) E 


96.744580 


12.188 


3.327x10" 


■6 


5 




2(1,1) —> 1(1,0) E 


96.755510 


20.108 


2.496x10" 


■6 


5 


"CO" 


1^0 


110.210353 


5.288 


6.389x10" 


■8 


3 


HNC" 


1 


90.663593 


4.350 


2.709x10" 


■5 


3 


HCNT 


li ^ Oi 


88.630416 


4.253 


2.444x10" 


-5 


3 




I2 ->0i 


88.631847 








5 




lo ^Oi 


88.633936 








1 


N2H+T 


li ^Oi 


93.171880 


4.471 


3.654x10" 


-5 


9 




I2 ^ Oi 


93.173700 








15 




lo ^ Oi 


93.176130 








3 


" Jk.Ki - 


> J^, p,, quantum numbers. 













^ J^J' quantum numbers. 

3p ^ Jp, quantum numbers. 
" Energies of the upper levels relative to the Jk = l-i level for E-type CH3OH, and 
relative to the Jk = Oo level for A-type CH3OH. 
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Table 2. Noise detection limit on the spectra with no emission. 



Source 


Tmb 3-0" detection limit (mK) 




CH3OH 


i^CO N2H+ HNC 


HON 


GO.32-0.20 


230 






G6.6 1-0.08 


104 






GIO. 10-0.72°= 


69 


414 142 127 


163 


G30.78+0.23 


95 







" Walsh et al. (1997) originally detected a 6.67 GHz methanol 
towards GIO. 10-0.72. Subsequent observations (S. EUingsen, private 
communication) have failed to detect the maser emission, hence we 
discard the source from ours sample in the remainder of this work. 



